NASA TM X- 



NASA TECHNICAL 

MEMORANDUM 


NASA TM X- 71962 

COPY NO. 



ANALYSIS OF THE FEEDBACK SYSTEM IN A NONINTRUSIVE 
DYNAMIC FUWEIER FOR l€ASURING F060 OSCIliATIONS 


/ioch 

WCEACK SYSTPM tu ANALYSIS OF THE 

CSCL 1 


By William G, Chapin 


N7«-2305tf 


63/ 1 4 


foiclas 

39385 


This informal documentation medium is used to provide accel erated or 
special release of technical information to selected users. The contents 
may not meet NASA formal editing and publication standards, may be re- 
vised, or may be incorporated in another publication. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
LANGLEY RESEARCH CENTER, HAMPTON, VIRGINIA 23665 



1. Report No. 2. Government Accession No. 

NASA TM X-71962 

3. Wecspier*i's Catalog No. 

4. Title and Subtitle 

ANALYSIS OP THE FEEDBACK SYSTEM IN A NONINTRUSIVE 
DYNAMIC FLOWMETER FOR MEASURING F0GO OSCILLATIONS 

5. Report Date 

6. Performing Organization Code 

7. Author (s) 

William G. Chapin 

3. Performing Organization Report No. 

10. Work Unit No. 

502-32-02-06 

9. Performing Organization Name and Address 

NASA, Langley Research Center 
Hampton, VA 23665 

11. Contract or Grant No. 

13. Type of Report and Period Covered 

Technical Memorandum 

12. Sponsoring Agency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546 

14. Sponsoring Agency Code 

15. Supplementary Notes 


Interim technical information release „ subject to possible revision and/or 
late formal publication 


16. Abstract 

Equations are developed which describe the closed loop feedback system 
operation of a proposed ultrasonic, dynamic, noiaiaitrusive flowmeter whose 
design is based on a constant phase, voltage controlled frequency feedback 
concept. These equations are based on linear feedback system theory. The 
time constant of s lew pass filter ic taken into account. The equations 
show that the larger the open loop gain, the smaller the error due to 
fluctuations in the speed of sound and tha mailer tho effective time 
cons tant. 


17. Key Words (Suggested by Auther(s)) (STAR category underlined) 

Instrumentation, flowmeter, 
nonintrusive, ultrasonic 


18 . Distribution Statement 


Unclassified - Unlimited 


Star Category 


19. Security Classif. (of this report) 

Unclassified 


20. Security Classif. (of this page) 
Unclassified 


21. No. of Paces 

43 


22. Price* 

$5.25 


* Available from i 


(The National Technical Information Service, Springfield, Virginia 22151 

(STIF/NASA Scientific and Technical Information Facility, P.0. Bok 33, College Park, MD 20740 

Q 

1 
















ANALYSIS OF THE FEEDBACK SYSTEM IN A NONB1TRUSIVE DYNAMIC FLOWMETER 
FOR MEASURING POGO OSCILLATIONS 


By William G. Chapin 
Langley Research Center 


SUMMARY 

Equations are developed which describe the closed loop feedback system 
operation of a proposed ultrasonic, dynamic, nonintrusive flowmeter whose 
design is based on a constant phase, voltage controlled, frequency feedback 
concept. These equations are based on linear feedback system theory. The 
time constant of a low pass filter is taken into account. The equations 
show that the larger the open loop gain, the smaller the error due to 
fluctuations in the speed of sound and the smaller the effective time 

constant . 


INTRODUCTION 

Many liquid-propellant rocket vehicles have experienced longitudinal 
vibrations because of an instability arising from interaction of the vehicle 
structure with the propulsion system. These vibrations, nicknamed Fogo 
after the jumping stick, have occurred principally in the first longitudinal 
structural mode during the first stage of a launch vehicle. The vibration 
begins spontaneously, intensifies, and then dies away - typically in a 
period of 10 to 40 seconds. 

In the interest of predicting whether or not the Space Shuttle Vehicle 
will "Pogo" it is necessary to measure the dynamic flow component in the 
LOX fuel lines during ground testing of the space shuttle main engine. This 
measurement should identify fluctuations in flow velocity to 100 Hz over a 
range +1.5 m/sec to +15 m/sec (±5 ft/sec to ^0 ft/sec) with a resolution 
of .003 m/sec (.01 ft/sec). 

gj_nce no flowmeter was available which coulo. meet the stringent 
requirements of this measurement, one had to be developed under an R&D design 
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study contract as reported in NASA CR-112313 (ref. l). This flowmeter 
operates on the principle of propagating two continuous ultrasonic sound 
waves (one upstream and one downstream) through the flowing medium by means 
of externally mounted transducer-receivers, (fig. l). This flowmeter 
apparently meets the requirements of nonintrusiveness , fast response, high 
accuracy and high resolution, and was designed to measure flow rates in 
cryogenic liquids and water. The heart of this flowmeter is its closed 
loop feedback system which represents the state-of-the-art in flow technology. 

Reference 1 contains an analysis of the flowmeter design concept, along 
with a block diagram, of the electronics and some circuit and mechanical 
details; however, reference 1 contains no analysis of the closed loop feed- 
back system operation. It was concluded that an analysis of the closed loop 
feedback system operation was needed to determine how or to what extent the 
system would actually produce the behavior described in the reference 1 
design analysis. Such an understanding was considered necessary to the 
effective evaluation of the design concept. For example, in reference 1, 
it is shown that the propagation of both an upstream wave and a downstream 
wave would, by means of a cancellation process, eliminate a measurement 
error due to fluctuations in the speed of sound in the flow medium. The 
question arose as to whether complete cancellation or only partial cancella- 
tion is an inherent feature of the actual systoa. It was realized that a 
mathematical model of the closed loop feedback system would be very helpful 
in answering this question. 

In this report, such a model is developed. It is based on classical 
linear feedback system theory and small fluctuations about a static 
operating point. 

First, the pertinent aspects of the analysis of reference 1 are 
presented. It shows that the simultaneous transmission of a downstream and 
an upstream wave results, ideally, in the elimination of fluid velocity 
measurement errors due to fluctuations in the speed of sound. Then the 
analysis of the closed loop servo operation is presented with appropriate 
servo block diagrams depicting the transfer functions and supporting 
equations. 
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SYMBOLS 


amplitude of the voltage controlled oscillator (VCO) output 


c 


velocity of sound 


D 


tube diameter 


f 


1 


y 


f 


10 



frequencies of downstream and upstream systems 9 respectively 

fixed components of the frequencies of the downstream and upstream 
systems, respectively 


F 


1 


9 


variable components of the frequencies of the downstream and 
upstream systems, respectively 


K & gain constant of dc amplifiers 


K, 


K 


m 


K l- K 2 


gain constant of low pass filters 

gain constant of multipliers 

gain constant of transmitter~receivers 

gain constants relating the voltage to the variable frequency 
component of the VCO units for the downstream and upstream 
systems, respectively 


K , K, , K , gain constants 

3 4 5 shown in figure 3. 


of the various flowmeter components as 


M, N integers 

ou t output at a block in figure 3, where n represents the number 

v of the block. 


S Laplace operator 

t time 

T time constant of low pass filters 


velocity of fluid 


V 


1 


9 


instantaneous outputs of the amplifiers of the downstream and 
upstream systems, respectively 


fl* 



instantaneous outputs of the low pass filters of the downstream 
and upstream systems, respectively 
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V t y instantaneous outputs of the multipliers of the downstream and 
101 m2 upstream systems, respectively 


-1 <4- ^ - -< 


^ U,,* 4- Ur 


v component of V . assumes not uO ue uOm^icucij atucuua^u w * AC 

^ low pass filter 

V component of V assumed not to be completely attenuated by the 

1112 low pass filter 

0 acute angle between the direction of wave propagation and the 

wall of the tube (see fig. l) 

<J> , ij> phase lags between the transmitter and receiver for downstream 

1 2 A oirotgmo respectively 

UiiVA. VtJMPW WA VJ w - y — --“f -- ^ 

X wave length of the transmitted ultrasonic wave for the downstream 

system 

subscript s static operating values about which systems operate dynamically 
subscript d. dynamic values about the sca&ic operating values 


superbar, ”, Laplace transform of the time domain dynamic component of a 
over a variable 
variable 
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GENERAL THEORY AND DESCRIPTION OF 
PROPOSED SYSTEM 


The basic purpose of the proposed nonintrusive dynamic flovmeter 
(ref. l) is to measure the dynamic flow velocity of liquids in a tube. In 
figure 1 each transducer contains a transmitter and a receiver. The trans- 
mitter of transducer 1 and the receiver of transducer 2 are called the 
’’downstream system.” The transmitter of transducer 2 and the receiver of 
transducer 1 are called the "upstream system.” Each of the transmitters 
is driven by a digital frequency synthesizer. Pertinent aspects 
of the analysis of the flowmeter design concept contained in reference 1 
are summarized next. The important outcome of. the analysis is that with 
suitable adjustment of some system parameters # the difference in voltage 
outputs of the two systems can be made a function of axial fluid velocity 
and independent of the speed of sound. This measurement concept represents 
the ideal behavior to be attained as closely as possible by the actual 
system. 

Consider a wave of frequency f propagated downstream from the 
transmitter of transducer 1 to the receiver of transducer 2 at an angle 9 
with the wall of the tube. The phase of the wave at 2 is 

2ttD . 1 
^1 sin 9 X 

also, 

c + v cos 9 

f i 

If these equations are combined, there is obtained 

2*1*1 (l) 

^1 = (c + v cos 9) sin 9 
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In the anticipated applications , it has been estimated that with a 

™3 

single frequency system 9 a change in pressure of only 10 atmospheres “would 

_ _ ___ * a ^ ^ a J im -4- ^ nwft ^-P PAim/1 onnVi + ivM cnpr 1 ! fn Pfitinn 

(jttUbC tx V^riablUJll X1J. uuw w V * ~jt — 

of .003 m/sec would be exceeded. Pressure variations are very likely to 
be considerably more than this, not to mention other sources of variation 
in the speed of sound. To minimize the effect of changing sound speed, a wave 
is also propagated upstream from the transmitter of transducer 2 to the 
receiver of transducer 1. Here 

2TTDf„ 

_ £ /p} 

*2 ” (c - v cos 0) sinO 

For the downstream system 9 an increase in fluid velocity tends to 
decrease the phase angle, while for the upstream system, such an increase 
tends to increase the phase angle. On the other hand, for each system, an 
increase in the speed of sound tends to decrease the phase angle. The 
frequency f^ consists of a fixed frequency f^ 0 plus a voltage controlled 
frequency The latter F-j^ is caused to vary so that 


(J) 1 = f (2N - 1) (3) 

This is an idealized relationship, since, as is shown later at least a 
small deviation of from ~ (2I~l) is an inherent feature of the servo 

system. 

Similarly, f consists of a fixed component, f 2Q plus a variable 
component Fg such that , ideally , 

<t> 2 = \ (2M - 1) W 

The frequency f is selected to be 1 MHs and fgQ is set at 1.1 MHz. 
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Equations (l) and (2) can be rewritten as 


27TD(f 10 + Y ± ) 

^1 “ (c + V cos 0) sin 0 


( 5 ) 


*2 (c - v cos 0) sin 0 


(6) 


If the right sides of (3) and (5) are equated and the resulting equation 


is solved for F^, there is obtained 


F, = iShil-SlaJL (. , v cos «) - f w 


(T) 


Similarly, from (4) and (6), there is obtained 


f 2 = ■ 3 - 4s JL - (0 - V cos e) - f 20 


( 8 ) 


By operating the digital frequency synthesizers in the frequency module' 
tion mode, F^ and F g can be made proportional to voltage inputs to the 
synthesizer . 

That is 


and 


F 1 * K 1 V 1 


F 2 K 2 V 2 


(9) 

(10) 


Furthermore, to obtain an output independent of the speed of sound, it 
is necessary to adjust and so that 
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( 11 ) 



Finally, from equations (?) through (ll) s an equation is evolved which relates 
a voltage output, — V^, to the axial fluid velocity vs 


V 


1 



(2N-1) 

*1 


sin 8 cos 6 
2D 


v 





( 12 ) 


Here the troublesome speed-of- sound variable cancels out and the last 
term is merely a dc bias. Equation (12) then represents the ideal behavior 
to be attained by the actual flowmeter system. 
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DEVELOPMENT OF THE LINEAR FEEDBACK MODEL 


Figure 2 is a simplified diagram of one channel of the proposed flowmeter 
(appendix G, ref. l). Notation is that pertaining to the downstream system 
which is analyzed in some detail. Figure 3 shows the more detailed block 
diagram of the flowmeter (ref. 1, page 57). The "Filter" of figure 2 
corresponds to the "LPF" (Low Pass Filter) of figure 3. An analysis is 
first made of the simplified system of figure 2, after which it is shown 
that the system of figure 3 essentially corresponds to that of figure 2. 

The relationships for the upstream system are easily obtained from the results 
of the analysis of the downstream system. 

On various occasions , use is made of the following trigonometric 

identities: 

sin a sin b = ^ cos (a - b) - ^ cos (a + b) (I.l) 

cos a cos b - — cos (a - b) + ^ cos (a + b) (1.2) 

2 2 

The 1/2 term does not explicitly appear in any equations since it can 
be included as part of a gain constant. 

Analysis of Downstream System 

For the downstream system (fig. 2), equations which express pertinent 
instantaneous time relationships are now restated or developed. 

The inputs to the multiplier are the reference wave and the reference 
wave delayed in phase by <]>^» and as (I.l) implies, 

V = K i [ A ain 2,r(F i + f io> *] [% A 8in ( 2 " (F i * V ‘ - *i)] 

- V* A 2 CO 5 -*! - K mV 2 005 + 2 f lo' 1 " *i] 

where K - K*/2 
m m 
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If it is assumed that, in passing through the filter, the high frequency- 
component is completely eliminated, the component of interest would become 

V = K K , A 2 cos <j>- (13) 

m^ m <p 1 

As shown previously, 

* „> 

T 1 sin 0 (c + v cos 0J 

It is planned to make the time constant, T, of the low pass filter 
about .01 seconds (ref. 1), which implies that the response is down 3 db 
at about 16 Hz. Therefore, the output of the filter is a function of the 
time variation of the input. That is. 



(t) = V 

*1 



(t)] 


(14) 


Since the maximum frequency of interest is 100 Hz, it is assumed that 
any time constants of the amplifier are insignificant. Hence, 

V = K V (15) 

1 “ f l 

and, as stated previously, 

F 1 = K 1 V 1 ^ 9 ) 

If small fluctuations in c and v about a static operating point 
are assumed, the dynamic relationships can be considered to be equivalent 
to the differential relationships. Then from equation (5) 
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and with 


F, = F. , v = v and c = c 
1 Is’ s s 


2ttD ... AV a- 27rD ■ ■ (V + f..Uc + v cos 6) (-l)dc 

d< frl = sin 0 (c + v cos IT 1 8in 6 ls 10 8 8 

s s 


+ fy + f )(c + v cos 0) (~l)(cos 0)dv 

sin 0 ls 10 s s 


d(t 1 = 


sin 6(c + v cos 6) 1 

s s 


2irD(F lg + f 10 ) 

11 “ ----- ^ 

sin 9 (c + v cos 6) £ 
s s 


2ttD(F + f, n ) COS e 

ls _10 dv 

2 

sin 0 (c + v cos 0) 
s s 


If the subscript "d" is used to denote the dynamic component, the 
preceding equation can be rewritten as, 


2ttD(F 1s + fpo) 


. . _ F /+\ — . c 

*ld ttJ = sin 0 (c g + v g cos 0) Id sin 0 + COB 0 ) 2 d 


3 S 


2ttD(F, + f, J cos 9 
iS 12 v d (t) 

sin 0 (c + v cos 0) 

5 S 


Likewise from equation (13) 


dV =-K K. A 2 (sin )d4>, 
Tn m 9 -1- 
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( 17 ) 


•which implies V^a^) = ” K m K <j, A ( sin 
Also from equations (14)* (15), and (9) 

V nd (t) = v fia [ \ia (t)1 

v ia (t) “ K a v na (t) 

F ia ct) " *1 v ia (t) 


(18) 

(19) 

( 20 ) 


Now take the Laplace transforms of equations (16) through (20) 


V s > - 


2ttD 


sin 6(c + v cos 0) 

8 S 


v s) 


2ifl)(F 


is 


f 10^ 


sin 0(c v cos 9) 
s s 


0 c.(S) 
2 d 


sin 9(c + v 

s s 


) cos 0 
cos 0)^ 


v d (S) 


(21a) 


or 


~^ld^ sin 


- 2ttD 

eTc + v cos 0 ) * Id. 
s s 


F. ,(S) + 


* V 


sin 0(c + v cos 0)‘ 

s s 


V s ) 




gTP(F i s * f io 
sin 0(c + v 

8 S 


) cos 0 
cos 0)^ 


v d (S) 


(21b) 


For the multiplier, the S domain equation for the dynamic voltage output 
is 


-v 

V 

mid 


(s> - v* A2(sin V (-?id ls)) 


( 22 ) 
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For the filter. 


K„ 


^fld^ = TS + 1 


Vd (s) 


(23) 


and for the output amplifier. 


*id ts) - K a (s) 


(2k) 


The S domain equation for the controlling frequency is 


f ia Cs) = K i ? id ls) 


(25) 


Equations (21b) through (25) can be represented by the block diagram 
shown in figure 4. It is desired to transform the configuration shown in 
figure 4 to a configuration of the form shown in figure 9* To do this, the 
general block diagram transformations shown in figure 5 are used. 

The diagram shown in figure 6 is derived from that shown in figure 4 by 
applying transformation T.l to the blocks in the forward loop and T.3 to the 
blocks imm ediately ahead of the inputs. The diagram shown in figure 7 is 
derived from that shown in figure 6 by applying transformation T.2 to the 
block to the left of the feedback loop. Finally, the block diagram shown in 
figure 8 is derived from that of figure 7 by use of transformation T.l. The 
diagram shown in figure 8 has the desired configuration. R in figure 9 
corresponds to (S) + (cos 0)v^(S) in figure 8. 

For a negative feedback system represented by the block diagram shown 
in figure 9, the relationship between the output C and the input R is 

C = r— ~ R, which implies that 
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K K K K A 2 (sin <t>. )2TTD(F n _ * f, J 

"a"r m <j> J-s 


V s) 


sin 0(c g + v g cos 0) (TS * 3.) 

LK K1K, A 2 (sin <j> n J2 ttD 
, . T. a f m 4> I s . 

1 + sin e(c + v cos 0}(TS -5- 1} 


(c^(s) + (cos 0)v^(S)) 


(26a) 


or 


Wm K * A2(sln V^^ls * f 10 : 


V s) * 


TS+1 + 


sin 0(c + v cos 

£ - b — — 

ILK K K K, A 2 ( sin (!) )2 itD 

1 a f m o is 

sin 0(c + v cos 0) 
s s 


(c d (S) + (cos 0)v d (S)) 


(26b) 


The question arises as to what is the sign of sin At present, it 

is assumed that the magnitude of the right hand term of the denominator of 
equation (26a) is to be made »1. It is shown later that this is a necessary 
condition for the desired operation of the flowmeter. Inspection of 
equation (26b) reveals that if sin 4> la is negative, there would be a pole in 
the right hand plane which would result in an unstable system. An analysis 
of this problem is presented next. 

It is assumed that no dynamic fluctuations are present but that, for some 

reason, $ is not at the static operating point If ^ were in either 

the second or the third quadrant, then cos ^ would be negative. Inspection 
of equations (13), (5), (lM, (15) and (9) reveals that the output of the 
multiplier would be negative, as would be V fl and, therefore, 

A negative V would cause to decrease which, in turn, would cause 
$ to decrease. Consequently, the system would drive toward the design 
Btetic operating condition ^ - (2» - 1) § . nhsre H is odd. Therefore, 
sin <(> would be positive. If, on the other hand, <J> X were in either the 
fourth or the first quadrant, then cos <j) x would be positive which would 
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cause F to increase. Consequently, <J> 1 would increase, and the system 
would again drive toward the condition 4^ = (2N - l) where K is odd. 
Therefore, sin 4> would again be positive. If, for some reason, 4>^ 
were to exactly equal (2N - l) where N is even, it seems likely that 
very soon some external or internal disturbance would occur and the 

nulling process would be triggered. 

In the remainder of this report, it is assumed that N and M (for the 

upstream system) are odd. 

The static operating point, 4>^ s * would not be exactly at (2N - l)— 
because a non zero cos 4^ is needed to produce at least a small multiplier 
output voltage, V^, so that the amplifier can maintain the required non 
zero variable frequency component F^ An equation is now developed which 
relates this difference between 4>^ s and (2N - l) — to system parameters. 
First one must rewrite equations (13) » (5), (1*0 » (15), and (9) in terms 
of values at the static operating point: 


V , = K K. A cos 4> 1e 

mis m <J> Is 


(27) 


4>. 


Is 


* V 

sin 0(c + v cos 0) 

s s 


(28) 


V fls V fls^mls^ K f V mls ’ 


since 


static flow conditions are assumed. (29) 


Furthermore 


V ls ' K a V fls 


(30) 


and 


F ls = *1 V ls 


(31) 


The difference between (2N - l) — and can be written 


= f (2» - 1) - 4> ls 


(32a) 
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which implies 


(32b) 


Is 2 


(235 - 1) - A<j> 


Then from equations (32a) and. (28) 


2uD F n 2 itD f ln 

7T ( v 13 j-U I'Z'l) 

A^>i = o’ (2N - 1) - s i n 0(c t v msT) sin 0(c * v cos 0) l J ' ' 

w s s s s 

BVrm. omio-H f\n ( . COS A = COS (~ (2S - l) “ A$ ) 02” COS <£-, „ = Sin A*, 

■ * " ^ is ’ £_! J- -LO - 1 

If sin is expanded in a Mac Lauren series 9 ‘there is obtained 




cos * la » Bio 4* x - 


Therefore, equation (27) becomes 


2 

Vs ' K m K * A (A h - IT 


From equations (3U), (29) » (30), end (31) 


2 ( ^1 } 

F is * ¥.W* A - “55 " • - •> 


The phase angle , <J) , attained by the system at ™ 0 is 


. = g » D f io 

P 10 sin 0(c * v COG 0) 

S S 


If equations (35) and (36) are substituted into (33) s it "becomeB 
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Comparison of the denominator on the right side of equation (37a) with 

that of (26a) reveals that it is 1+ the open loop gain constant, except for 

the sin di term of the latter equation. Assume, more specifically. 

Is 

that the system is designed such that 


2it D K, K K,. 
1 a f 

sin 0 (c + v 

S 


K K 


m 

cos 


s 


£ 



> > 1 


which implies that the second term on the left hand side 
of equation (37a) 



Therefore 




At most , f (2H-1) - <I> 10 = ± * radians. 


max 


i + 


TT 

2ttD K.K K,K K, A 2 
1 a f at 9 

sin 0(c + v cos 0) 

S 3 


(37b) 


Further, the right hand side of equation (3Tb) is <<1, yielding the 
relationship 


( A ^) 3 

3! 


<< A<j> 1 


and equation (37b) can be written 


(A$ ) = — « 1 (37c) 

1 I max 27TD K.K K-K K, A 

j + 1 a f m (p 

sin 0(c + v cos 0) 

s s 
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which implies that 


♦lb - 


(38a) 


sin 4> ls = 1 


(38h) 


Equation (38a) verifies that, for static conditions, the idealized 
condition stated by equation (3) will be closely attained. 

Simplified expressions for equations (26a) and (26b) can be obtained; 
From equations (28) and (38a) 


= 


2 ” D(r is ♦ V 

sin 0(c + v cos 0) 

s s 


£ I (2N - 1) 


which implies 


( F l g + f 1Q ) 2 ^ 2N L 


sin 6(c + v cos 0) 

s s 


2ttD 


(2N - l) sin 0(c + v cos 0) 

S 3 

!+D 


(39) 


The block diagram obtained by applying equations (38b) and (39) to 
figure 8 is shown in figure 10. 


Equations (26a) and (26b), which relate system output to inputs, now become 

. 2 , 


K a K f K m K 4> A (2N " 1)Tr 
2(c + v cos 0)(TS + l) 


V s > ■ 


2TO K.K KJC K. A 
1 a r m q> 

^ + sin 0{c + v cos 0j"(TS + 1) 

s s 


(c d (S) + (cos 6)v d (S)) 


(l*0a) 


or 
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w/ (2N - 11 l 

2(c + v cos 0 ) 

_ s s /_ 

v ld (s) = ( c d^ SJ + ( cos 6 > v d (S) ) 

2ttD K,K KJC K, A 2 
TS + 1 + — 1 a " f m 

sin 0(c + v cos 0) (bob) 

s s 


The denominator is of the form 

IT 1 

t me* 

TTY” + x \ 


This shows that the effective time constant of the closed loop 
system is reduced by the factor 1 ■> K, where K is the open loop gain constant. 

Analysis of Upstream System 

For the upstream system, the analysis makes use of the instantaneous 
equations : 


V = K K , 
m2 m <j) 


*2 " 


f2 


COB <j> 2 


(hi) 

ro 

o 


(6) 

- V cos 

0) 

( v (t) 

) 

(U2) 


V 2 = K a Y f2 


F 2 = K 2 V 2 


(hi) 

( 10 ) 
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The analysis of this upstream system is similar to that of the 
downstream system j however, whereas in the downstream system, an increase m 
fluid velocity tends to cause a decrease in phase angle, in this upstream 
system, an increase in fluid velocity tends to cause an increase in phase 
angle. Thus, 


21TD F 2d (t) 

4u d U) ■ " 

d sin 0(c - v cos 0) 

s s 


2 " D(F 2S + v 


sin 0(c - v cos 0) 

s s 


2 ' c d (t) 


2ttD (F 0 + f^) (cos 0) 

+ 2= 22 . v d (t) 

sin 0(c - v cos 0) 

s s 


In the S domain 


= - 


2ttD F 2d (S) 
sin 0(c - v cos 9) 

5 S 


g * D < F 2 e + W 

sin 0(c - v cos 0) 

s s 


2 • ‘d (S > 


2TTD(F gs - f 2 „)(d °s 8) _ 

— ■ ~ 2 ~ > 

sin 0(c - v cos 0) 

s s 


The other Lap lace- transformed dynamic equations are: 


W s) 

* V/ sin *2s (-+2d (S>> 

(46) 

W s) 

- Kf • W S) 

TS + 1 

(47) 

V 2d (S) . 

‘ Vf2d (S) 

(48) 

F 2d (S) ■ 

■ K 2 V 2d (S) 

(49) 


Equations (45) through (49) can he put into the block form shown in 
figure 11. 
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Figure 11, "by transformations similar to those used in the downstream 


system , assumes the tuim shewn in .figure 12. 


Also, by analogy with equations (38b) and (39)» the following relationships 
are used: 


F + f 
2s 20 


sin $ 2s = 1 
sin 0 


(2M “ l)(c “V cos 0) 
4D ss 


Then, from the general output/input relationship for a negative feedback system, 

wa — : - 1 (= a ( s > - 0 > ^ (s) > 

2(0 - V cos e) (TS + 1) 

V ls) - S s 3 < 50 > 


2d 


2ttD K„K KJC K, A 
2 a f m <p 

+ sin 0(c - v cos 6)(TS + 1) 

s s 


Analysis of Combined Downstream and Upstream Systems 

The block diagrams for the two individual systems of figures 10 and 
12 can be combined as shown in figure 13. In the feedback block of the 
upstream system, has been expressed in terms of by use of equation (ll). 


From equations (UOa) and (50) 


V s) - V s) = 


K KJC K. A (2H - 1 )JI _ _ 

2 (c H cos e)(ns * i) (c d (s) * <cos e)T d (s) > 

s a 

2ttD ILK K^K K. A 2 
1 a f ml 

^ + sin 0"(c " + v cos 0)(TS ^T) 
s s 

K KJC K, A 2 (2M -1)tt _ _ 

zCc ' - v' col WtSwT" ( =d (S) " (oOS 0)v d (S)) 


1 + 


2tiD KjyC^K, A‘ 

sin 0(c - v cos 0)(TS + l) 

s s 


( 51 ) 
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Next, it is shown that, if the open loop gain (GH in figure 9) is 
> > 1, and, if and Kg are adjusted for the relationship expressed by 
equation (ll), the ideal output/input relationship expressed by equation (12) 
is closely attained. If equation (12) which expresses an instantaneous 
relationship is rewritten to express a dynamic relationship, there is obtained. 


V id< fc > - 


V 2d (t) = 


- /2N_zi\ 

'1 kJ 


sin 9 cos 8 
2D 


v d (t) 


(52) 


Equation (52), then, expresses the relationship between dynamic 
YoXtage output and dynamic fluid velocity that would be obtained by an ideal 
system. Equation (12) and therefore, equation (52) assume that and Kg 
are related in the manner expressed by equation (ll). 

jf -kjjg right hand terms in each of the denominator expressions are 

> > 1 , 

» ld (S) - V 2d (S) i p ° ln '■ ( 'd (S) + (C0S 61 ”d (S)> 


. (gM 4P 1 * Bln 6 ' <c d <s> - (cos 0) 7 d (s)> 


If equation (ll) is used for Kg 


V. .(S) - V (S) = - 1 ) sin g- (c (S) + (cos e)v (s) ) 

ld Zd kD ^ 


or 


. I 21 * rir ( ”d (s) ’ !o ° b e) 7 d (s,) 

V - V fs) = ~ 1 ) sin 6 c< ^ s ^ v (S) 

V, j(S) v 2d^ b ‘' ?n k d 


Id' 


2D K., 


(53) 


In the time domain, 


V ld^ “ V 2d^ 


(2N - 1) sin 9 cos 9 y 


2D K. 


(5^) 
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As seen by inspection of figure 13 and equation (5^), the system out- 
put is inversely proportional to the feedback path gain, 4DK^/(2N - l)sin 0. 
This, of course, is an inherent disadvantage in the use of negative feedback. 
To obtain large open loop gain, and at the same time, adequate output/input 
sensitivity, adequate forward path gain is needed. Changes in the static 

levels of fluid velocity and sound speed, v and c , change the forward path 

s s 

gain, but high open loop gain reduces this effect. 


ANALYSIS OF THE BLOCK DIAGRAM OF THE PROPOSED FLOWMETER 

Next, it is shown that the proposed system of figure 3 corresponds 
essentially to the simplified system (fig. 2) upon which our analysis 
has been based. 

For the sake of simplicity, interactions between the upstream receiver 
and downstream transmitter and between the downstream receiver and upstream 
transmitter are neglected. 

In figure 3, the signal outputs are shown at various key points in 
the system. Some of these output equations are developed below where it is 
felt they might not be easily determined by inspection. 

The outputs are denoted by cut/^N where "n" represents the number 
of the component block shown in figure 3. 


°ut0= k; ( K 3 a sin 271.^1 ) (K l ,K 2 A sin [27Tf 2 t - <fr 2 - a]) 

f 

By use of trigometric identity (I.l) and the relationship, K^= K^/2 

Out(g) = A 2 [cos (2-rrf 2 t - <J > 2 - a - 27Tf ^t ) 

- cos ( 2iTf 2 t - <f> 2 - a + 2irf^t ) ] 


2k 



or 

o “ t ®' 
Similarly 


K K^ A 2 |cos(2ir[f 2 - f 1 ]t - <J> 2 - a) -cos(2Tr[f 1 + f 2 ]t - <t> 2 - a)} 


°utg^= A 2 [cos(2TTf 2 t - a - Sirf^t + -cos(2uf 2 t - a + Znt^ -4^)] 


or 


Out fix- a 2 |cos(2Tr[f 2 - f-Jt - a + 4^) -cos(2Tr[f 2 + f-Jt - a - 4^)} 


and Out /^\ is Kg x the ( f 2 - f^) term of Out (l^ 


Similarly 


Out 


1? ,= KgK^K^ A 2 cos(2-rr[f 2 - f-Jt - a) 


Out 


0 = K^KgK^K^ A 2 cos(2TT[f 2 - t ± ]t -4> 2 -a^KgK^ A 2 cos(2ir[f 2 -f-^t-a)} 


i . 

By use of trigometric identity (1.2) and the relationship, = K^/2, 
Out(g)= KgK^K^ k k |cos(2Tr[f 2 - - <t> 2 - a + 27r[f 2 - - a) 

+ cos(2ir[f 2 - f x ]t - a - 2ir[f 2 - f^t + a + <t>^ 


or 


Out0= KgK^K^ A^ |cos(2tt . 2[f g - fjt - <P 2 - 2ot) + cos 4> 2 J. 


Similarly 


Out0= A^ | cos(2tt * 2[f 2 - t ± ]t - 2a + 4^) + cos 
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It can be seen from the foregoing equations that, except for constant 

terms: 

1. Equation (9) holds for both the actual system and the simplified system 

2. The unfiltered part of the output of ^2) corresponds to equation 

(13). 

3. The output of (23) corresponds to equation (14). 

4. The output of (£5) corresponds to equation (15). 

5. Equation (5), of course, holds for both the actual and simplified 
systems . 

Similarly, for the upstream system, equation (10) holds for both the 
actual system and the simplified system and : 

1. The unfiltered part of the output of (2^ corresponds to equation (4l). 

2. The output of (25) corresponds to equation (42). 

3. The output of (26) corresponds to equation (43). 

4. Equation (6), of course, holds for both the actual and the simplified 
systems. 


CONCLUDING REM&3KS 


An analysis of the feedback system used in servoing the transmitting 
frequencies of the Panametrics, Inc. dynamic ultrasonic flowmeter has been 
made. The analysis is based on the assumption that small fluctuations in the 
speed of sound and axial fluid velocity occur about a static operating point 
representing a steady-state fluid flow rate. 

The analysis reveals that the design concept is essentially sound as 
long as the open loop gain can be made »1. This condition is needed to 
cancel out errors due to fluctuations in the speed of sound and to decrease 
the flowmeter's effective time constant needed to obtain that frequency 
response to at least 100 Hz. For example, aa open loop gain of 100 results 
in an effective time constant of 10 ^ seconds. Further, adequate forward 
gain is needed to obtain the desired output/imput sensitivity, and at the 
same time, obtain the high open loop gain required for stability. Unf or- 
nately, the forward gain varies with variations in the static (steady-state 
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flow rate) operating values of fluid velocity and sound speed, v g and c g , 
respectively, but high open loop gain reduces this effect also. 

Based on the root locus stability criteria, the flowmeter's eervoing 
was found to contain only a single pole in the left hand plane and is 
inherently stable. 

As a result of this analysis, the flowmeter's innovative concept, it is 
felt, will be more clearly understood and its ability to measure oscillatory 
flows more effectively evaluated. 


27 



REFERENCES 


1. Pedersen, Norman E. , Lynnworth, Lawrence C., and Carnevale, Edmund H. : 

Nonintrusive Dynamic Flowmeter. Final Report, June 19T3, NAS1-11T56 , 
NASA CR - 112313 

2. Ogata, Katsuhiko: Modern Control Engineering. Prentice-Hall, Inc., 

1970, pp 92-93. 


28 






vco 

# 1 


A sin 2n(F, + f in )t 



Figure 2. Simplified operational diagram of flowmeter. 
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Figure 6* Block diagram derived from figure 
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Figure 7 . Block diagram derived from figure 6 . 
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Figure 8 


Downstream system block diagram derived from figure 7. 
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Downstream system block diagram derived from figure 8 by use of 

relationships sin 0. =1 and F, + f ln = (2N - l)sin 9(c + v cos 9)/4D* 

is is 1U s s 


Figure 10. 
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Figure 11. Block diagram representation of S domain 
equations for upstream system. 
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Figure 12. Upstream block diagram derived from figure 11. 












